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Abstract 

We present a search for the rare flavour-changing neutral-current decay B~ —* K~vV based on 
a sample of (86.9 ± 1.0) x 10 6 T(45) — > BB events collected in the BABAR experiment at the 
SLAC B-factory. Signal candidate events are selected by fully reconstructing a B + — > D°X + 
decay, where X + represents a combination of up to three charged pions or kaons and up to two 
7T° candidates. The charged tracks and calorimeter clusters not used in the B + reconstruction 
are required to be compatible with a B _ — ► K~vv decay. We observe a total of three signal 
candidate events with an expected background of 2.7 ± 0.8, resulting in a preliminary limit of 
B(B~ — > K~vT>) < 1.05 x 10 -4 at the 90% confidence level. This search is combined with the 
results of a previous and statistically independent preliminary BABAR search for B~ — > K~vV to 
give a limit of B(B~ — > K~vv) < 7.0 x 10 -5 at the 90% confidence level. 
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1 INTRODUCTION 



The quark level process b — > svv represents a rare flavour-changing neutral-current (FCNC) de- 
cay which proceeds at the one-loop level in the Standard Model (SM) via "penguin" and "box" 
diagrams such as those shown in Fig. ^ The inclusive b — » svv process is nearly free from theoret- 
ical uncertainties associated with strong interaction effects, permitting a fairly precise prediction 
of the SM branching fraction. The inclusive branching fraction, summed over the three neutrino 
flavours, is estimated to be (4.1 jli'o) x 10~ 5 _F. Since additional heavy particles would also con- 
tribute additional loop diagrams, various "New Physics" scenarios can potentially lead to significant 
enhancements to the SM branching fraction Unfortunately, an experimental search for the in- 
clusive b — > svv process is extremely difficult in a B-factory environment due to the presence of 
two unobserved neutrinos which limit the available kinematic constraints that can be exploited in 
order to suppress other B decay backgrounds. 




Figure 1: Electroweak penguin (left) and box (right) Feynman diagrams for the process b — > svv 
predicted by the SM. In both cases the amplitudes are expected to be dominated by the heavy t 
quark contribution. 

Instead, we search for the exclusive B~ — ► K~vv decay mode, which proceeds via the b —>■ svv 
process. The SM branching fraction for B~~ — > K~vv is estimated to be B(B~ — > K~vV) ~ 4 x 
10~ 6 [3J|I]. The best published limit on the exclusive branching fraction is from CLEO [5] with a 
limit of B(B~ —> K~vv) < 2.4 x 10~ 4 at the 90% confidence level. BABAR has already reported a 
preliminary upper limit B(B~ -> K~W) < 9.4 x 1(T 5 6 based on 50.7 fb" 1 data. The two BABAR 
analyses use reconstruction methods which produce mutually exclusive data samples, permitting 
the two statistically independent results to be combined to obtain an improved limit. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring 
during the period from 2000 - 2002, corresponding to a total integrated luminosity of 80.06 fb _1 
collected at the T(4S) resonance. This "onpeak" sample is estimated to contain (86.9 ± 1.0) x 10 6 
BB pairs. This data set is supplemented by a sample of 9.6 fb _1 of "offpeak" data collected 
approximately 40 MeV below the T(4S) resonance, which is used to study continuum background 
sources due to e + e~ — > // where / = u,d, s, c, e, fi, r. 

The BABAR detector is an hermetic detector optimized to provide precision vertexing, charged 
and neutral particle reconstruction and particle identification in an asymmetric B-factory environ- 
ment. Tracking is provided by a five-layer double-sided silicon vertex tracker (SVT), surrounded 
by a 40-layer drift chamber (DCH) filled with a mixture of helium and isobutane. The SVT and 
DCH are situated within a 1.5 T solenoidal field. K/n separation is provided by a quartz ring- 
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imaging Cherenkov detector (DIRC) located immediately outside of the DCH. The electromagnetic 
calorimeter (EMC) is used to measure energy and position of photons and electrons. Muon iden- 
tification is achieved through segmentation and instrumentation of the iron of the magnetic flux 
return (IFR) using resistive plate chambers. A more detailed description of the BABAR detector 
can be found in |7j. 

A GEANT4 [H] based Monte Carlo (MC) simulation is used to model the signal efficiency and 
physics backgrounds. MC samples equivalent to approximately three times the data luminosity were 
used to model BB events, and samples equivalent to approximately 1.5 times the data luminosity 
were used to model continuum events. 

3 ANALYSIS METHOD 

Due to the presence of two unobserved neutrinos, the B~ — > K~uu decay mode 7 lacks the kinematic 
constraints which are usually exploited in B decay searches in order to reject both continuum and 
BB backgrounds. Consequently, the approach which is used in this analysis is to first reconstruct 
the accompanying "tag" B + , which is produced in association with the signal B~ through the 
process T(4S) — > B + B~ , and then search for evidence of a B~~ — > K~vu decay among the tracks 
and clusters not associated with the reconstructed tag B. In order to avoid experimenter bias, the 
signal region in data is not examined ("blinded") until the cuts are finalized. 

We reconstruct the tag B + in a set of decay modes B + — > TP X + where X + is a hadronic 
system composed of up to three charged mesons (either ir or K) and up to two 7r° candidates. The 
D candidate is reconstructed in one of the three decay modes D —* K + -k~, D — » K + tt~it or 
~ir + ir . Candidate B + decays are identified by combining D candidates with sets 
of charged tracks and n° candidates until the combination yields a B candidate consistent with 
the kinematics expected for a true B meson decay. We use the two kinematic variables tties and 

AE, defined by = \J -Ebeam 2 — Pb 2 where pb is the momentum vector of the B candidate and 
-E'beam is the beam energy, and AE = Eb — -Ebeam) where Eb is the energy of the B candidate. 
All these quantities are evaluated in the center of mass (CM) frame. If multiple B candidates are 
identified within the kinematic region m-Es > 5.2 GeV/c 2 and —1.8 GeV < AE < 0.6 GeV, only 
the candidate for which AE is closest to zero is retained. 

Combinatorial backgrounds from continuum processes are significantly reduced by requiring the 
thrust, computed using all tracks and clusters in the event, be less than 0.925 and that |cos#t|) 
the magnitude of the cosine of the angle between the thrust axis defined by tracks and clusters 
used to reconstruct the tag B candidate and the thrust axis defined by all other tracks and clus- 
ters in the event, be less than 0.8. Correctly reconstructed B meson candidates produce a peak 
in the mES distribution above a combinatorial background at the nominal B mass as shown in 
Fig. [2j The tag B yield is determined directly from data by determining the peaking component 
of the mES distribution. Tag B candidates which are reconstructed in MC events containing a 
true B~ — > K~vv decay are found to possess very little combinatoric background, since there are 
few additional tracks and clusters in the event which can be randomly combined to produce com- 
binatoric tag B candidates. Events with tag B candidates lying within a signal region defined by 
5.272 GeV/c 2 < toes < 5.288 GeV/c 2 are retained for use in the search for B~ — > K~vv decays. 
Events in the region 5.225 GeV/c 2 < toes < 5.265 GeV/c 2 are retained for use in background 

7 Charge conjugate modes are implied throughout this paper, however the signal mode will always be denoted as 
a B~ decay, while the fully reconstructed tag B will be denoted as a B + decay to avoid confusion. 
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Figure 2: The distribution m-Es for B + — » D X + candidates in the B~ — » K~vv signal MC (top), 
and in the data (bottom). The bottom plot also shows the expected contributions from continuum 
and BB MC. No signal-side selection cuts have been applied. Signal MC is shown scaled to the 
data luminosity assuming B{B~ — > K~vV) = 4 x 10~ 6 ; inclusive background MC and offpeak data 
are shown scaled to the onpeak data luminosity of 80.06 fb _1 . 



studies as discussed below. A discrepancy of approximately 25% between the observed yield in 
data and the predicted yield in MC is corrected by scaling the peaking component of the MC 
simulation. Good agreement between data and MC is obtained once this correction is applied. The 
same scale factor, (75 ±7)%, is used to correct the tag B yield in signal MC as well as MC estimates 
of peaking backgrounds. The quoted error of 7% is due to the uncertainty in the estimation of the 
peaking component of the data and due to the discrepancy between MC and data in the shape of 
the non-peaking background component. 

Once a reconstructed tag B candidate has been identified, B~ — > K~vV signal candidates are 
selected by considering all tracks and clusters in the event which are not used in the tag B recon- 
struction. This set of tracks and clusters is referred to in the following as the "signal-side" of the 
event. Candidate events are required to possess exactly one signal-side reconstructed charged track 
with a charge which is opposite that of the tag B. The signal-side charged track multiplicity is 
plotted for signal MC and data in Fig. |3l The signal candidate track is required to satisfy parti- 
cle identification criteria for charged kaons based on information from the tracking system and the 
DIRC. The kaon candidate is boosted into the CM frame assuming a kaon mass hypothesis, and the 
CM momentum, p* K , is required to be greater than 1.5 GeV/c. The average particle identification 
efficiency in the momentum range of interest is ~ 85% and the typical tt — > K misidentification rate 
is ~ 2%. The p* K of signal candidate tracks is plotted in Fig. 0] We assume the kaon momentum 
spectrum is described by the decay model of reference and correct the signal MC distribution 
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Figure 3: The distribution of the number, Arracks, of signal-side charged tracks is plotted for the 
B~ — > K~vV signal MC (top), and for onpeak data and generic MC (bottom) for events which pass 
the tag B selection. No signal-side selection cuts have been applied. Signal MC is shown scaled to 
the data luminosity assuming B{B- -> K~vV) = 4 x 1(T 6 , and generic MC and offpeak data are 
shown scaled to the onpeak data luminosity of 80.06 fb _1 . 



(which is generated with a phase-space model) accordingly. 

Constraints are also imposed on signal-side EMC clusters to reject events with significant neutral 
energy deposition in the calorimeter. B~ — » K~vv events possess an average of two additional 
signal-side clusters. These are generally attributable to hadronic split-offs in the EMC, usually 
from pions or kaons associated with the tag B side of the event, or to beam related backgrounds. 
Events possessing one or more 7r° candidates, composed of combinations of two EMC clusters 
with lab frame energy greater than 30 MeV which combine to produce an invariant mass in the 
range 122 MeV/c 2 < m 77 < 145 MeV/c 2 , are rejected. In addition to this 7T° veto, a restriction is 
imposed on the total "extra" signal-side neutral energy, E e ^x x&) that is present in the event. -Eextra 
is computed by summing the CM-frame energies of all signal-side EMC clusters with lab frame 
energy greater than 30 MeV. Signal candidate events are required to possess E'extra < 300 MeV. 
The .Eextra distribution obtained from signal MC is compared to data in Fig. 03 

The requirement that signal candidate events possess low signal-side track multiplicity and 
relatively little additional neutral energy in the EMC also tends to select B + B~ events in which the 
tag B + has been correctly reconstructed, but the opposing B~ has one or more unreconstructed 
particles which have passed outside the detector acceptance in either the forward or backward 
directions. An additional cut is therefore imposed on the direction of the missing momentum 
vector in the CM-frame, p^iss) which is required to satisfy | cos^p^J < 0.8. 

Due to the low signal-side multiplicity, there is almost no tag B reconstruction mode-dependence 
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Figure 4: The CM momentum distribution of signal-side kaons candidates is plotted for the 
B~ — > K~vV signal MC (top), and for onpeak data and generic MC (bottom). Plotted events are 
required to possess exactly one signal-side track satisfying kaon identification criteria and having 
a charge opposite that of the tag B. Signal MC is shown scaled to the data luminosity assuming 
B(B~ — > K~vV) = 4 x 10~ 6 , and generic MC and offpeak data are shown scaled to the onpeak 
data luminosity of 80.06 fb _1 . 



of the measured signal-side efficiency. The overall selection efficiency, etot , can therefore be expressed 
as the product of the tag B selection efficiency, etag, and the signal-side efficiency e s i g . The overall 
efficiency is estimated to be e to t = (0.046 ± 0.005)%, while the e s i g is estimated to be (35 ± 1)%. 
The uncertainties on the efficiencies are due to both statistics and systematics. 

Monte Carlo modeling of the signal efficiency and background estimates is validated by com- 
paring the yields obtained in a number of data control samples. These samples include 9.6 fb _1 of 
offpeak data, an toes sideband region defined by 5.225 GeV/c 2 < toes < 5.265 GeV/c 2 , an -Eextra 
sideband defined by 0.5 GeV < E extTll < 1.5 GeV, and a "large toes" sideband spanning the entire 
region defined by 5.225 GeV/c 2 < toes < 5.265 GeV/c 2 and -Eextra < 1-5 GeV. In addition, we 
retain samples of events which pass all of the nominal signal selection requirements except that 
they are required to have a total of two or three charged tracks associated with the signal side of 
the event instead of only one. All data control samples were found to be in good agreement with 
the MC predictions, as shown in Tabled 

Backgrounds consist primarily of B + B~ events in which the tag B + has been correctly re- 
constructed but in which the accompanying B~ decays to a high- momentum kaon and additional 
particles which are not reconstructed by the tracking detectors or calorimeter. Typically these 
events contain one or more K\ and/or neutrinos, and frequently also one or more additional 
charged or neutral particles which pass outside of the tracking or calorimeter acceptance. This 
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Figure 5: The total extra neutral energy distribution, -E ex tra is plotted for the B~ — > K~vv signal 
MC (top), and for onpeak data and generic MC (bottom). Plotted events are required to possess ex- 
actly one signal-side track satisfying kaon identification criteria and having a charge opposite that of 
the tag B. Signal MC is shown scaled to the data luminosity assuming B{B~ — > K~vT>) = 4 x 10 -6 , 
and generic MC and offpeak data are shown scaled to the onpeak data luminosity of 80.06 fb _1 . 



Table 1: A comparison of data yields and MC predictions in the signal region and in various control 
samples. Quoted uncertainties reflect MC statistics only. 



MC type 


Signal Region 


"lES 


large m E s 


-^extra 


^tracks — 2 


-^tracks — 3 


B+B- 
B°W 
uu, dd, ss 
cc 
T + T - 


1.7 ±0.6 






1.1 ±0.5 


1.8 ± 1.0 
1.8 ± 1.0 



7.0 ± 1.4 

I. 4 ±0.6 
14.0 ±2.9 

II. 1 ±2.6 




3.3 ±0.9 
0.6 ±0.4 

2.4 ± 1.2 
2.4 ± 1.2 




17.4 ± 1.9 
0.9 ±0.5 
0.6 ±0.6 
1.9 ± 1.0 



54.6 ± 3.4 
3.5 ± 1.0 
1.2 ±0.9 
3.1 ± 1.4 



Onpeak data 
Total MC 
(£ = 80.06 ftr 1 ) 


3 

1.7 ±0.6 


7 

4.8 ± 1.7 


31 

33.5 ±4.2 


10 

8.8 ±2.0 


21 

20.7 ±2.3 


55 

62.4 ±3.9 


Offpeak data 
Continuum MC 
(£ = 9.6 fb- 1 ) 




0.11 ±0.05 




0.4 ±0.2 


1 

3.0 ±0.5 




0.6 ±0.2 




0.3 ±0.1 


1 

0.5 ±0.2 
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"peaking" background component is evaluated directly from B + B~ MC and estimated to yield 
1.7 ± 0.6 events in 80.06 fb _1 of data. This estimate is validated by comparison with the data in 
the -Eextraj -^tracks = 2 and A^ ra cks = 3 sideband, all of which have large contributions from peaking 
B + B~ backgrounds. Due to the limited MC statistics, a smaller combinatorial component of the 
background is estimated by scaling the observed yield in the MC tties sideband into the signal re- 
gion. This scaling assumes an "Argus function" shape for the tties distribution of the combinatorial 
component which is obtained from data using a "wrong-sign" tag B sample in which the charge 
of the D° daughter kaon is inconsistent with the charge of the tag B candidate, resulting in no 
significant peaking component in the tties distribution. Scaling the toes sideband into the signal 
region yields an additional background of 1.0 ±0.4 events, leading to an estimated total background 
of 2.7 ± 0.7, where the quoted uncertainties are due to MC statistics. 

4 SYSTEMATIC STUDIES 

Estimates of systematic uncertainties are summarized in Table [2j Systematic uncertainties in the 
branching ratio determination are dominated by the statistical uncertainty on the background 
estimate, and by the uncertainty on the determination of the tag B yield in data using (75 ±7)% 
scaling correction determined from the comparison of B + D X + yields in data and MC. bmce 
the tag B yield is determined directly from data, no other systematic uncertainties associated with 
the tag B need to be assigned. The signal track reconstruction and kaon identification procedure 
yield comparatively small uncertainties on the efficiency and background estimates. Uncertainties 
also arise from the MC modeling of the energy and multiplicity of low energy clusters in the EMC 
which would potentially produce a bias in the -Eextra distribution. These uncertainties are estimated 
by evaluating the change in the MC efficiency and background estimates when low energy clusters 
are selectively removed from the MC until data and MC multiplicity distributions are in agreement. 
An additional systematic uncertainty is assigned to the efficiency estimate as a result of the re- 
weighting procedure used to correct the MC kaon momentum spectrum to be consistent with that 
predicted by reference [3]. Note however that this uncertainty does not account for variations in 
the momentum spectrum which would result from the use of other theoretical models. A small 
uncertainty also enters the branching ratio limit calculation from the estimation of the number of 
B + B~ events present in the data sample. 

5 PHYSICS RESULTS 

Unblinding the analysis revealed a total of AT se i = 3 events in the signal region, with an expected 
background, A^, g , of 2.7 ± 0.8 where the additional systematic uncertainties from table|2]have been 
combined with the MC statistical uncertainty. The mES distribution and signal kaon candidate 
momentum spectrum are shown in Fig. Eland Fig. respectively. The B~ — > K~vv branching ratio 
is computed as follows: 

B(B- _ K-uV) = (1) 
z ■ 1\ B + B - • etot 

where e to t = (0.046 ± 0.005)% is the overall signal selection efficiency, and N B + B ~ is the number 
of T(4S) —* B + B~ events in the data. N B + B - = (43.5 ± 0.5) x 10 6 is obtained by assuming equal 
branching fractions for T(45) decays into charged and neutral B mesons. The central value of 
the branching ratio is determined to be B(B~ — > K~vV) = (0.8 ± 2.0) x 10 -5 where the quoted 
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Table 2: Systematic uncertainties 



Source 


Relative uncertainty 


Signal efficiency: 




Signal MC statistics e s i g 


1% 


Signal MC statistics et ag 


5% 


Tag B yield 


7% 


Signal track reconstruction efficiency 


1% 


K particle ID efficiency 


2% 


-^extra 


2% 


Kaon momentum correction 


3% 


Total cre/e 


10% 


Background estimate: 




Generic MC statistics 


27% 


Tag B correction 


7% 


Track efficiency 


5% 


-E"extra 


8% 


Total background estimate uncertainty 


29% 


B + B- yield 


1.2% 



uncertainty is from systematics only. The significance of the central value is somewhat less than la 
and we therefore quote a 90% confidence level limit. The branching ratio limit is computed using 
a frequentist approach based on reference 0- The confidence level for a given branching fraction 
limit "guess" is obtained by generating a large number of experiments in which the systematic 
uncertainties in the inputs to equation ^ are modeled by Gaussian distributions and the signal 
statistics are modeled by a Poisson distribution. The limit is set at the value of the branching 
fraction at which 10% of the generated experiments produce a yield which is less than the observed 
data yield of three events. This procedure results in a limit of B(B~ -> K~uV) < 1.05 x 10" 4 at 
the 90% confidence level assuming the central model of reference !3_. This limit is weakly model 
dependent, since the signal efficiency depends on the kaon momentum spectrum. Varying the model 
within the range specified in [3j results in a variation in the limit in the range (1.02 — 1.10) x 10 -4 . 
Alternatively, the model of reference fl] yields a limit of 9.5 x 10 -5 . Sensitivity to new physics 
will also be model dependent, however it should be noted that this analysis has an efficiency which 
is relatively uniform as a function of kaon momentum above the 1.5 GeV/c cut. The efficiency is 
zero for kaon momenta below 1.5 GeV/c, which may effect the limit interpretation for exotic New 
Physics modes with significantly different kaon momentum spectra. 

This analysis can be combined with the result of a preliminary BABAR search for B~~ — » K~vv 
reported previously [5], based on a sample of 50.7 fb _1 of data. This analysis yielded a limit of 
B(B~ -> K~uV) < 9.4 x 10~ 5 at the 90% confidence level assuming the model of reference [2]. 
In this analysis, the tag B was reconstructed in a set of semileptonic B decay modes of the form 
were X can be either nothing, or one or more photons which are consistent with 
the decay products of higher mass open charm states such as D*° — > D°^/tt°. Two events were 
observed, and were treated as signal for the limit determination (i.e. no background subtraction 
was performed). Because this previous analysis required an identified high momentum lepton for 
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Figure 6: The tues distribution of events passing all other B 
three selected events in the signal region. 



K vv selection cuts, showing the 



> 

<D 



© 
O 



■4-" 

c 

> 

LU 




Onpeak data 
Total MC 
B*B" MC 
B°B° MC 
cc MC 

uu, dd, ss MC 

x*x MC 



2.5 



Signal track momentum (GeV/c) 



Figure 7: The CM momentum distribution of signal kaon candidates which pass all other 
B~ — > K~vV selection cuts, showing the three data events above the 1.5 GeV/c cut. 



the tag B reconstruction, while in the present work the tag B modes are purely hadronic, the 
two analyses are by construction statistically independent and can be readily combined. Using the 
same frequentist method as described above and ignoring correlated systematic uncertainties which 
are expected to be small, we obtain a combined limit of B(B~ — > K~vv) < 7.0 x 10~ 5 at the 90% 
confidence level. 



6 SUMMARY 

We have performed a search for the rare FCNC decay B~ — > K~vv using a method in which the 
accompanying tag B meson is reconstructed into a set of hadronic final states. Using a data sample 
corresponding to an integrated luminosity of 80.06 fb -1 , we observe a total of three signal candidate 
events, consistent with the background expectation of 2.7 ± 0.8. We determine a preliminary limit 
of the branching fraction B(B~ — > K~vv) < 1.05 x 10 -4 at the 90% confidence level. Combining 
our result with the result from a previous and independent BABAR search for B~ —* K~vv yields 
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a combined preliminary limit of B(B — ► i/z/) < 7.0 x 10 5 . The result is consistent with the 
SM expectation of B{B~ K~vv) ~ 4 X 10~ 6 . 
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We present a search for the rare flavour-changing neutral-current de- 
cay B~ -> K~vV based on a sample of (86.9 ± 1.0) x 10 6 T(45) -» 
events collected by the BABAR experiment at the SLAC 5-factory. Signal 
candidate events are selected by fully reconstructing a B + — > D X + de- 
cay, where X + represents a combination of up to three charged pions or 
kaons and up to two 7r° candidates. The charged tracks and calorimeter 
clusters not used in the B + reconstruction are required to be compatible 
with a B~~ — > K~W decay. We observe a total of three signal candidate 
events with an expected background of 2.7 ± 0.8, resulting in a prelimi- 
nary limit of B(B~ — > K~vv) < 1.05 x 10~ 4 at the 90% confidence level. 
This search is combined with the results of a previous and statistically in- 
dependent preliminary BABAR search for B~ — > K~vV to give a limit of 
B(B~ JTw) < 7.0 x 10~ 5 at the 90% confidence level. 
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